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Abstract 


We propose a soluble model to incorporate the nonlinear effects in the tran¬ 
sition probabilities of the multiphonon Giant Dipole Resonances based on 
the SU(1,1) algebra. Analytical expressions for the multi-phonon transition 
probabilities are derived. Enhancement of the Double Giant Resonance exci¬ 
tation probabilities in relativistic ion collisions scales as {2k + l){2k)~^ for the 
degree of nonlinearity {2k)~^ and is able to reach values 1.5 — 2 compatible 
with experimental data. The enhancement factor is found to decrease with 
increasing bombarding energy. 
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Coulomb Excitation in collisions of relativistic ions is one of the most promising methods 
in modern nuclear physics M- One of the most interesting applications of this method 
to studies of nuclear structure is the possibility to observe and study the multi phonon 
Giant Resonances [|I|. In particular, the double Dipole Giant Resonances (DGDR) have 
been observed in a number of nuclei The “bulk properties” of the one- and two- 

phonon GDR are now partly understood and they are in a reasonable agreement with 
the theoretical picture based on the concept of GDR-phonons as almost harmonic quantized 
vibrations. 

Despite that, there is a persisting discrepancy between the theory and the data, observed 
in various experiments |]6HlO| that still remains to be understood: the double GDR excitation 
cross sections are found enhanced by factor 1.3 — 2 with respect to the predictions of the 
harmonic phonon picture |^, 0, [^, [|^]. This discrepancy, which almost disappears at high 
bombarding energy, has attracted much attention in current literature [Q, i2^-p4|; 

among the approaches to resolve the problem are the higher order perturbation theory 
treatment [^, and studies of anharmonic/nonlinear aspects of GDR dynamics |^, p2|,^, 
m- Recently, the concept of hot phonons [|^, [|^] within Brink-Axel mechanism was 
proposed that provides microscopic explanation of the effect. These seemingly orthogonal 
explanations deserve clarihcation which we try to supply here. 

The purpose of this work is to examine, within a a soluble model the role of the nonlinear 
effects on the transition amplitudes that connect the multiphonon states in a heavy-ion 
Goulomb excitation process. Most studies of anharmonic corrections [p2|-|2^ concentrated 


on their effect in the spectrum . Within our model, the nonlinear effects are described 


by a single parameter, and the model contains the harmonic model as its limiting case when 
the nonlinearity goes to zero. We obtain analytical expressions for the probabilities of 
excitation of multiphonon states which substitutes the Poisson formula of the harmonic 
phonon theory. For the reasonable values of the nonlinearity, the present model is able to 
describe the observed enhancement of the double GDR cross sections quoted above. 

Having in mind to show how analytical results follows from the nonlinear model, and 
to explain how the model works, we restrict ourselves here to its simplest version (trans¬ 
verse approximation, or SU(1,1) dynamics) and keep numerics up to minimum level. We 
postpone till further publications detailed numerical analysis and comparison with the data. 
Microscopic origins of the nonlinear effects (which are considered here phenomenologically) 
are also beyond the scope of this presentation. 

We work in a semiclassical approach []n| to the coupled-channels problem, i.e., the 
projectile-target relative motion is approximated by a classical trajectory and the excitation 


of the Giant Resonances is treated quantum mechanically 


The use of this method 


is justihed due to the small wavelenghts associated with the relative motion in relativistic 
heavy ion collisions. The separation coordinate is treated as a classical time dependent 
variable, and the projectile motion is assumed to be a straight line |^. 

The intrinsic dynamics of excited nucleus is governed by a time dependent quantum 

n 


Hamiltonian (see Refs. 


). The intrinsic state |'0(t) > of excited nucleus is the 


solution of the time dependent Schrodinger equation 

= |ffo + r(«))l IV'(*) > (1) 

where Hq is the intrinsic Hamiltonian and V is the channel-coupling interaction. We use 
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the system of units where h = 1, c = 1. The standard coupled-channel problem for the 
amplitudes a„(f) reads 

>= ^n{t) \n > exp ( - iEntj, (2) 

where En is the energy of the state \n > in the wave packet In our treatment, the 
nuclear states are specihed by the numbers of excited GDR phonons, N or n. Taking scalar 
product with the states < A^|, we get the set of coupled equations for the amplitudes a„ as 
functions of impact parameter b 

ia^{t) = J2< N\V\n > (3) 

n=0 

We assume the colliding nuclei to be in their ground states before the collision. The am¬ 
plitudes obey the initial condition ^ —cxo) = 6{n, 0) and they tend to constant values 
as t —>■ ±cxo (the interaction V{t) dies out at t — >■ ±cx)). The excitation probability of an 
intrinsic state |A^ > in a collision with impact parameter b is given as 

WN{b) = |ajv(oo)|l (4) 

and the total cross section for excitation of the state |iV > is given by the integral over the 
imact parameter 


ctat = 27r J bWN{b)db (5) 

bgr 

with the grazing value bgr = 1.2(Rg/^ -|- A\/p) as the lower limit. Hereafter, the labels exc 
{sp) refer to the excited (spectator) partner in a colliding projectile-target pair. We neglect 
the here nuclear contribution to the excitation process. 

In the following, it is convenient to treat the coupled channel equations ([^),(^ in terms 
of the unitary operator f// (in the interaction representation) that acts in the reference basis 
of multiphonon states (including the ground state |0)): 

4t/,(t) = Vi{t)Uj{t), Vi{t) = Ui{t = -oo) = I, (6) 

at 

and the time-dependent Hamiltonian H{t) = HQ + V{t) that acts in the intrinsic multi-GDR 
states is given by 


Ho = UjNd, Nd = Y^ dl^ndm, 

m 

V{t) = Vi{t)[{El_i)'^ - (El+iY] + vo{t){Eloy + Herm.Conj. (7) 

where El\^ and Elm are the dimensionless operators acting in the internal space of the multi- 
GDR states. In the harmonic approximation, they are given by the GDR phonon creation 
and destruction operators of corresponding angular momentum projection m, El\^ = d]^. 
The function Vi{t) is given by 
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Vi{t) = 


w 


[1 + (^t)2]3/2’ 


W = p 


ZspC^'y 


^excZexc 


2&2 \ AlLlmN-S^MeV' 


( 8 ) 


(the corresponding expression for no(t) can be fonnd in Ref. |^). Here, mj^ and e are the 
proton mass and charge, Z, N and A denote the nuclear charge, the neutron number and 
the mass number of the colliding partners, 7 = (1 — is relativistic factor, v is the 

velocity and the parameter p is the deal of the strength absorbed by the collective motion 
(usually assumed to be close to unity) [^. 

The harmonic approximation (ideal bosons) yields the transition probabilities between 
the states with numbers of GDR phonons differing by unity to grow linearly oc N. This 
model of ideal bosons has well known exact solution (see, e.g. [T^) that is given by the 
Poisson formula for the excitation probabilities 


I ^harm 12 


Wn = e- 


\^harm 12 


^harm 12 ^" 

W\ 


E l harm \ 2 
\^m \ 

m=0,±l 


r\\^harm\2 i \^harm\2 

I + I^Q I 


( 9 ) 


where the amplitudes are expressed in terms of the modihed Bessel functions. At 


\a 


harm 12 


0 


IS sup- 
“transverse 


the colliding energies sufficiently high, the longitudinal contribution (oc 
pressed by a factor proportional to 7 “^ |p . In the following, we will work in the 
approximation” dropping the longitudinal term for the sake of simplicity. This is a good 
approximation at high energies and the results are still qualitatively valid at lower energies. 

Our idea is to keep the spectrum of GDR system harmonic with the Hamiltonian Hq = 
ujN. That is supported by the systematics of the observed DGDR energies, E 2 , which yields 
E 2 — (1.75 — 2)u |l|. The conclusion on the weak anharmonicity in the spectrum follows 
also from theoretical considerations 


25, m 


The transition operators El\El can however include nonlinear effects. (In particular, 
this could be a result of nonlinearities in the phonon Hamiltonian obtained in higher orders 
of perturbation theory). A reasonable accounting for these nonlinear effects that we adopt 
in this work consists of generalization of the transition operators 


Z^ln — El 


2k 





( 10 ) 


where the parameter k ^ 0 determines the strength of the nonlinear effects, i.e., the problem 
reduces to the harmonic oscillator with linear coupling when 


1 


0 . 


( 11 ) 


If A; < 0, the transition probabilities are suppressed as N grows. Positive values of k that 
we will consider here correspond to the enhancement of the matrix elements. 

It is convenient to introduce the three following operators D^, D~ and 




1 


(d+i 


dt,)^j2k + N,, 

1 


= - 


(dX 


+1 


^ ~ 27/2 ~ '^- 1)1 

d~^i){d^i — d-i) + 2{2k + N()) , 


( 12 ) 
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with Nd = d^id^i + d1iid-i. It is easy to check that they obey the standard commutation 
relation for the noncompact SU(1,1) algebra 


=£)-, L)+,L)° =-£>+, D-,D- 


= 2D°. 


(13) 


The dynamics of the our system, in transverse approximation, can be written in terms of 
the operators , D~ and (^) only. In the interaction representation, the evolution 
equation (j^) takes the form 




^/k 




Vk 




(14) 


for Ui{t) (the last equation follows from (^), (0) after using the commutation relations (P!3[ ) 
and [Nd,D^] = ±D^. 

From purely mathematical viewpoint, the problem described by the last equation drops 
into the universality class of the systems with SU(1,1) dynamics that can be analyzed by 
means of generalized coherent states f^for the SU(1,1) algebra. In particular, the problem 


of the parametric excitation of a one-dimensional harmonic oscillator ( [^|, |^|) belongs to 
the same class. The physical meaning of the algebra generators (|^ in our case is of course 
quite different from that of p^j30[| . (For other algeabraic approaches to scattering problems, 
see Refs. |^, [0). 


The formal solution of our problem is given by the expression for the unitary operator 
Ui{t) as a time-ordered exponential (see, e.g., |^) 


Ui(t) = Texp 



(15) 


Due to closure of the commutation relations between the operators D^,D~ and that 
enter the exponential in Eq.([T5|), the time-ordered exponential can be represented in an¬ 
other equivalent form that involve ordinary operator exponentials only (see, e.g., |^). The 
operator Ui{t) can be expressed as 


Ui{t) = exp 


a{t) 

y/k 


D- 


exp 


logil- 


k 


- #(f) 




exp 


a*{t) 

y/k 


D- 


(16) 


where one has the product of ordinary exponentials which contain some time-dependent 
complex number a{t) (star means complex conjugation) and real number 0(t) (phase). The 
functions a{t) and 0(t) can be found from simple differential equations which relate the 
unknown a{t) and 0(f) with the function ni(f) in the Hamiltonian H{t). These equations 
(see below) can be restored after substituting the right hand side of Eq. (pT]) into the left 
hand side of the Schrodinger equation for the operator Ui{t) (0 and collecting the terms 
which have the same operator structure. 


Proceeding this way, we obtain, after some algebraic manipulations, from ([T^) with using 
the commutation relations (0), the formula 
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for differentiating the operator exponential and the Baker-Hausdorf relations for the operator 
exponentials [0, the following Riccati-type equation for the complex amplitude a: 

d c? 

, « = + . (17) 

dt k 

The expression for the phase (j){t) is given by a simple integral 0(t) = 

t 

(2//c) / dtiRe{vi{ti)a{ti)e~^‘^^^) (as the phase does not contribute to the squared abso- 

— OO 

lute values of amplitudes that enter the excitation probabilities and the cross sections, we 
will not be interested in it in the following). 

Once the solution for the differential equation ([T7| ) is found, the expression for the am¬ 
plitudes ajv(t) which we are interesting in follows from (pTf) immediately after projection of 
the state 


mt)) = r/,(«)| 0 > 


onto the intrinsic states with dehnite number of GDR phonons, N. For example, we have 
From Eq.(^) while acting onto the ground state a particularly simple result 


Ui{t)\0)=e 




k 


exp 


a{t) 

y/k 


D- 


| 0 ). 


We thus obtain the expression for the amplitude and the probabilities of the transitions from 
the ground state to the excited states with N phonons: 


|a(cxo)| 


= j ( j 


F(2A: +A^)V^^ /|a(cx))|2\^/^ 


Wn = |aAr(oo)p 


(18) 


Here, the quantity a(+cx)) is the asymptotic solution to the Riccati equation (|^ at t > cx) 
subject to the initial condition a(—cx)) = 0. The simple nonlinear equation ([T7|) abbreviates 
all orders of quantum perturbation theory for the problem Eqs. ( 1 . 0 ,( 0 ). It is also seen 
from Eq. ( 0 ) that unitarity is automatically preserved within present formalism {Uj = Uj^). 

The harmonic limit of these results corresponds to the case k ^ oo, when the nonlinearity 
disappears in the transition operators ( 0 ) and the coupling to electromagnetic held becomes 
linear (harmonic approximation). Then at fc —> cxo the last nonlinear term drops from 
the equation (P!7|), and the solution to the equation |q((cx))| reduces to its harmonic value 


\a 


harm \ 
±1 


given by modihed Bessel functions [0, [|I|, At the same time, the expression for the 


probabilities W, Eq. (piBf) , reduces at /c —> cx) to the Poisson formula (^, thus the harmonic 
results 0, i are restored. 

The simple nonlinear equation (0) for the amplitude a{t) can be readily solved nu¬ 
merically. One can construct the solution by means of successive approximations, provided 
that the nonlinearity is small enough to keep the hrst terms in the expansion to powers of 
\vi\/Vk, 


T Z 

a{t) = y dti + i J 


k 


iujti 


dti 


Vl{t2 




dU 


(19) 
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In fact, the first term of this expansion is an excellent approximation in many cases, especially 
for light projectiles/targets when |ni| (H) is small. We will use this approximation below as 
a good reference point to analyze the enhancement factor for excitation of DGDR. Keeping 
thus the hrst term in ([T^) we have 



00 

a(cx)) ~ 

—i J vi{t)e'^'^^dt 

— 00 


= 2-em) 

(j 


( 20 ) 


i.e., the well known expression expression in terms of the modihed Bessel function. While the 
amplitude a is thereby reduced to the harmonic value in this approximation, the excitation 
probabilities Wv are still given by the formula for the nonlinear model, Eq. Using 

(p!8|) and (pOD, one has therefore the approximate expression for the excitation probabilities 


^i( 0 |VV"r( 2 fc + iV) f2\K,{0?s 


2 „ 2 ' 


Wn = [I- 


N 


k 


mT{2k) 


2k 


( 21 ) 


here, s = 2—f '^ and f = —. The cross sections for the excitation of Wphonon GDR are 


vy 


(jjv — 27r 


T{2k + N) 
N\T{2k){2k)^ 


l^i( 0 |V' 

k 


dh 1 


2k 


mmvs 


2 „ 2 \Ar 


( 22 ) 


At nonzero nonlinearity l/{2k) > 0, the excitation probabilities for Multiple GDR 
turn out to be enhanced as compared to their values in the harmonic limit. In 


general, deviation in the excitation probabilities for the iV-phonon states Wn (p!8|) from 
their harmonic values are given by the ratio 



k 2k + i'] 

|a(+oo)|2^2fc 

a(-|-cx))p'^ 

TT/Ziarm J. X 1 

i=0 

1 y 

1 

|Q,fearm(_^QQ)|2iV 


In perturbative regime, when the amplitudes a (p!9|) are close to the harmonic solutions, the 
two last factors in this equation are close to unity, and the deviations are basically due to 


N-l 


the factor H is convenient to introduce the following ratio R 2 that characterizes 


i=0 


enhancement of the excitation probability of the two-phonon GDR due to nonlinear effects 
with respect to the transition probability in the harmonic limit 


R 2 = 


W 2 / Wi 




/ 


(23) 


From equations 


,(PkD we obtain the result 

2k + 1 |q;(-|-cx))P 


R 2 = 


2k 




(24) 


The hrst factor in this expression, results from the kinematic enhancement of the tran¬ 
sition probabilities due to nonlinear effects considered here. When the amplitude |q!(-|-cx))| 
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is small, the second and the third factors in Eq.(|^) are close to unity and we obtain the 
result 


2fc + 1 


2k 




(25) 


The cross sections can be obtained from the usual formula (^. For the ratio of the cross 
sections of excitation of the double GDR calculated within present nonlinear model ct 2 and 
within the harmonic phonon model we obtain in the same scaling regime 


r2 


0-2 


^harm 

^2 


V’ 


2fc + 1 
2k 


(26) 


where the exact numerical factor that stems from equations (PD, (|TD, (H), is close to 
unity when |q;| ~ -Cl. It contains dependence on the bombarding energy and 

describes deviations from the scaling behavior which are valuable at both low and high 
energies (we discuss these deviations below). 

The interesting feature of these results is that they are rather insensitive to many of 
the parameters that describe the colliding system. In the scaling regime, the enhancement 
factor r 2 depends on the parameter 1/k only, which is assumed to be an intrinsic property 
of the excited partner. As a consequence, the enhancement must be roughly independent 
of the properties of the spectator partner, once we study excitation of double GDR in the 
same nucleus. 

This is just what has been observed in experiments: the values of found for DGDR 
in 208 projectile using different targets U 0 are very close to each 

other and they correlate, within the error bars, with the value 


r 2 (^°®P 6 ) ~ 1.33 


(27) 


According to the scaling prediction, this corresponds to the nonlinearity parameter k equal 
to 


kr^pb) ~ 1.5 


(28) 


Numerical evaluation shows that the value of -0 (see Eq.(^)) is close to 0.88, thus, the value 
of the required nonlinearity is a bit bigger, corresponding to 1 /fc ~ 1 . 

The same picture was found in experiments on Coulomb desintegration of Au target 
using various projectiles Bi pj. Also, the similar conclusion of nearly 

constant value of r 2 in ‘^^^Pb target with various projectiles has been made in work fTH 


(though for differing experimental set-up). 

Below, we present the exact results for the cross sections calculated according to Eqs.(||), 
and (|I8| ) and with solving Eq.([T7|) numerically. 

The dependence of the enhancement factor r 2 = (T 2 / ([26|) for the DGDR excitation 
on the strength of the nonlinearity k is shown on Fig. 1. for the process 208 p^_j_ 208 p^ ^gg 
the case of bombarding energy e = 0.64GeR/per nucleon). It is seen that the enhancement 
factor drops to unity at big values of k (harmonic limit) and grows at stronger nonlinearity. 
The scaling value of r 2 is also shown for comparison. 

It is interesting to trace energy dependence of the enhancement factor. 







For typical projectile-target pairs and for moderate bombarding energies (7 ~ 1-3), both 


a 


and 


\a 


harm \ 
1 


are small as compared to unity and |a| ~ \a 


harm \ 
1 


The value of 


IS 


well approximated by and this results in approximate scaling behavior for the cross 

sections: one has r 2 = ct 2 /ct 2 “'’”* — For reasonable values of fc, the enhancement factor 
r 2 of the cross-section of Double Giant Resonance excitation can reach values ~ 1.5 that is 
compatible with experimental data |l[]. 

Deviations from this simple scaling rule occur at both low and high energies. At 7 —> 1, 


adiabatic approximation is valid, and this yields to |q:| > 


l^harm] 


Thus, R 2 > R 


scaling _ 


{2k -|- l)/{2k). By contrast, at higher energies, the dynamical nonlinear effects tend to 


reduce the magnitude of |a|, thus |q:|/|q: 


harm I 

1 I 


< 1 , and Ro < 


^scaling 


To sum up, the 


enhancement factor for the DGDR excitation cross section, r 2 = a 2 /<J 2 ^ drops from 2 — 
2.5 (for low bombarding energies £ ~ lOOMeR per nucleon) to 1.2 — 1.3 (for £ ~ 640 — 
700MeR per nucleon) while hxed value of k is used. On Fig.2., we plotted the value of 
the enhancement factor calculated numerically for the case of Pb process. The 

magnitude of nonlinearity is kept hxed, k = 1.3. 

These results are is in correspondence with the experimentally observed trends and with 
microscopic models based on Axel-Brink concept . 

To conclude, we presented here a simple model that accounts for the nonlinear effects 
in the transition probabilities for the excitation of multi-phonon Giant Dipole Resonances 
in Goulomb excitation via relativistic heavy ion collisions. The model is based on the 
group theoretical properties of the boson operators. It allows to construct the solution 
for the dynamics of the multi-phonon excitation within coupled-channel approach in terms 
of the generalized coherent states of the corresponding algebras. The well known exactly 
solvable harmonic phonon model appears to be a limiting case of the present model when 
the nonlinearity parameter 1/k goes to zero. The main advantages of the limiting harmonic 
case (unrestricted multiphonon basis, preservation of unitarity and possibility of analytical 
treatment) remain present in our nonlinear scheme. Therefore, the model can be viewed 
as a natural extension of the harmonic phonon model to include the nonlinear effects in a 
consistent way while keeping the model solvable. 

At high enough projectile energies, the problem is simplihed and the dynamics is governed 
by the generators of SU(1,1) group. The problem of the excitation amplitudes is reduced to 
single Riccati-type equation for the complex amplitude a. Its solutions can be obtained by 
means of successive approximations. In most cases relevant to applications, the perturbation 
is weak enough to be considered in the hrst order; the explicit solution for the amplitude is 
then given by the modihed Bessel functions. The probabilities to excite double-phonon GDR 
appear to be enhanced by an approximately universal factor {2k + l)/( 2 /c); the same scaling 
is roughly valid for the cross sections. This can be viewed as a hint that the discrepancy 
between the measured cross-sections of double GDR and the harmonic phonon calculations 
can be resolved within present nonlinear model by means of using an appropriate value of 
the nonlinear parameter 1/k for a given nucleus. The experimental values of enhancement of 
(72 with respect to the harmonic results for the excited nucleus are almost insensitive 

to the details of the collision process. On the other hand, the enhancement factor drops as 
the bombarding energy grows. Though the present results are obtained in the “transverse 
approximation” (neglecting the longitudinal response), they approximately hold in general 
case. This is consistent with the data and gives results similar to those recently obtained in 
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a possibly different context, with a theory based on the concept of fluctuations (damping) 
and the Brink-Axel mechanism |^, []^|, |^|. It would be certainly worthwhile to establish 
possible connections between the two approaches. 

The work has been supported by FAPESP (Fundacao de Amparo a Pesquisa do Estado 
de Sao Paulo). 
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Figure Captions 


Fig.l. Enhancement factor r 2 = for the Double GDR excitation in 208 p^_l_ 208 p^ 

process at bombarding energy e = 640Mel//per nucleon as a function of the parameter k 
(symbols, solid curve is to guide the eye). Dashed curve corresponds to the scaling r 2 ~ 

Fig.2. Enhancement factor r 2 = for the Double GDR excitation in 208 p^_l_ 208 p^ 

process as a function of relativistic factor 7 (symbols, solid curve is to guide the eye). The 
value of the nonlinear parameter is kept to be equal to k = 1.3. The scaling value (constant 
^ 1 ^) is shown by dashed curve. 
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